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Abstract 


Presented  in  this  report  is  an  in-situ  method  to  calibrate  X-band  antenna  arrays  desipied  to 
detect  and  track  low-level  missiles  over  water.  These  arrays  are  indispensable  in  ship  defence, 
because  they  can  operate  in  rainy  and  foggy  conditions  where  atmospheric  absorption  is  too 
high  for  the  use  of  more  accurate  millimetre  wave  radars  or  optical  devices. 

The  method  uses  a  set  of  array  snapshots  to  construct  a  calibration  matrix  for  the  antenna 
array.  It  incorporates  the  compensations  for  unequal  gain  and  phase  responses  at  the  array 
elements,  antenna  tilt,  mutual  coupling  among  the  elements,  the  effect  of  local  structures,  and 
small  errors  in  element  positions. 

Simulation  and  experimental  evaluations  showed  that  the  calibration  (a)  improved  the 
accuracy  of  target  elevation  angle  estimates  obtained  with  a  high-resolution  direction-finding 
method,  (b)  improved  the  resolution  of  target  and  image  elevation  angles,  (c)  suppressed  a 
bias  in  elevation  angle  estimates,  and  (d)  suppressed  some  anomalous  tracks  between  the  true 
target  and  image  tracks. 

The  calibration  improves  a  ship's  defence  against  low-level  missiles. 


Resume 

Ce  rapport  presente  une  methode  d’etalonnage  sur  place  de  reseaux  d’antennes  a  bande  X 
utilisees  pour  la  detection  et  la  poursuite  de  missiles  volant  a  basse  altitude  au-dessus  de 
I’ocean.  Ces  reseaux  sont  indispensables  pour  la  defense  des  navires,  car  ils  peuvent  operer 
dans  des  conditions  de  pluie  et  de  brouillard  pour  lesquelles  I’attenuation  atmosphdrique  est 
trop  grande  pour  I’utilisation  de  radar  millimetrique  plus  precis  ou  d’equipement  optique. 

La  methode  utilise  un  ensemble  de  mesures  instantanees  du  reseau  pour  construire  une 
matrice  d’etalonnage  pour  le  reseau  d’antennes.  Elle  incorpore  une  compensation  pour  les 
reponses  inegalees  en  phases  et  gain  au  niveau  des  elements  du  reseau,  de  I  inclinaison  de 
I’antenne,  du  couplage  mutuelle  entre  les  elements,  pour  les  effets  des  structures  locales,  et 
des  petites  erreurs  dans  les  positions  des  elements. 

Les  evaluations  experimentales  et  simulees  montrent  que  I’etalonnage  (a)  ameliore  la 
precision  des  mesures  de  Tangle  de  site  des  cibles  obtenu  avec  une  m6thode  de 
radiogoniomdtrie  de  haute  resolution,  (b)  ameliore  la  resolution  des  angles  de  site  de  la  cible 
et  de  son  image,  (c)  elimine  une  erreur  sur  la  mesure  des  angles  de  site,  et  (d)  elimine 
quelques  anomalies  entre  les  vraies  pistes  des  cibles  et  leurs  images. 

L’etalonnage  amdliore  la  defense  des  navires  contre  les  missiles  volant  a  basses  altitudes 
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Executive  summary 

This  report  describes  the  in-situ  calibration  of  X-band  antenna  arrays  designed  to  detect  and 
track  low-level  sea  skimming  missiles.  X-band  antenna  arrays  are  indispensable  in  ship 
defence  against  these  missiles,  especially  in  foggy  and  rainy  conditions  where  optical  and 
millimetre  wave  radars  become  ineffective  in  detecting  these  targets. 

The  antenna  arrays  must  be  calibrated  accurately,  because  a  more  accurately  calibrated  array 
yields  a  more  accurate  estimate  of  missile  elevation.  Practical  considerations  require  the 
calibration  be  carried  out  in-situ  with  a  remote  pilot  source.  This  report  describes  the  first 
successful  attempt  to  use  this  type  of  source  in  calibration. 

Simulation  and  experimental  evaluations  showed  that  the  calibration  (a)  improved  the 
accuracy  of  target  elevation  angle  estimates  obtained  with  a  high-resolution  direction-finding 
method,  (b)  improved  the  resolution  of  target  and  image  elevation  angles,  (c)  suppressed  a 
bias  in  elevation  angle  estimates,  and  (d)  suppressed  some  anomalous  tracks  between  the  true 
target  and  image  tracks. 

The  calibration  method  can  be  used  to  calibrate  the  ship-bome  X-band  Active  Phased  Array 
Radar  (APAR)  being  jointly  developed  by  Canada,  Germany  and  the  Netherlands. 
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Sommaire 

Ce  rapport  decris  I’etalonnage  sur  place  de  reseaux  d’antennes  a  bande  X  utilisee  pour  la 
detection  et  la  poursuite  de  missiles  volant  a  basse  altitude  au-dessus  de  I’ocean.  Les  reseaux 
d’antennes  a  bande  X  sont  indispensables  a  la  defense  des  navires  contre  ces  missiles, 
specialement  dans  des  conditions  de  pluie  et  de  brouillard  ou  les  instruments  optique  et  les 
radar  millim^triques  sont  inefficaces  pour  detecter  ces  cibles. 

Le  rdseau  d’antennes  doit  dtre  precisement  etalonner,  car  un  reseau  6talonner  produit  une 
mesure  plus  precise  de  I’angle  de  site  des  cibles  Des  considerations  pratique  exigent  que 
r^talonnage  soit  fait  sur  place  avec  une  source  pilote  a  distance.  Ce  rapport  decris  le  premier 
essai  reussi  pour  1’ utilisation  de  ce  type  de  source  pour  I’etalonnage. 

Les  evaluations  experimentales  et  simulees  montrent  que  I’etalonnage  (a)  ameliore  la 
precision  des  mesures  de  I’angle  de  site  des  cibles  obtenu  avec  une  methode  de 
radiogoniometrie  de  haute  resolution,  (b)  ameliore  la  resolution  des  angles  de  site  de  la  cible 
et  de  son  image,  (c)  elimine  une  erreur  sur  la  mesure  des  angles  de  site,  et  (d)  elimine 
quelques  anomalies  entre  les  vraies  pistes  des  cibles  et  leurs  images. 

La  methode  d’etalonnage  peut  etre  utilisee  pour  etalonner  les  radars  navals  a  antenne  reseau 
active  a  bande  X  developpee  conjointement  par  le  Canada,  L’Allemagne,  et  les  Pays-Bas. 
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1 .  Introduction 


X-band  antennas  are  indispensable  in  ship  defence  against  low-elevation  missiles. 
These  antennas  can  operate  in  rainy  and  foggy  conditions  where  atmospheric  absorption  is  too 
high  for  the  use  of  more  accurate  optical  devices  or  millimetre  wave  radars. 

Low-angle  tracking  (LAT)  methods  for  X-band  antenna  arrays  must  deal  with  the 
presence  of  interfering  image  signals  produced  by  reflections  off  the  water  surface.  They  can 
be  divided  into  two  groups.  The  first  group  uses  established  high-resolution  direction-finding 
(HRDF)  methods  such  as  the  MUSIC  method  [1]  and  the  Tufts-Kumaresan  singular  value 
decomposition  (TKSVD)  method  [2,3]  to  estimate  target  elevation  angles.  An  example  of 
application  is  given  by  Hayldn,  Greenlay  and  Litva  in  [4].  The  second  group  comprises 
several  maximum  likelihood  estimation  (MLE)  methods.  It  includes  (i)  the  Correlation  Height 
Analysis  (CHA)  method  and  its  variants  [5-10]  developed  by  Litva,  Rook,  Lo  and  Wong  and 
(ii)  the  refined  maximum  likelihood  (RML)  method  developed  by  Bosse,  Turner,  Lecours, 
Brookes  and  Riseborough  [11,13].  These  MLE  methods  use  a  highly  deterministic  model  of 
the  propagation  medium  to  model  the  interference  pattern  of  the  target  and  image  signals  at 
the  array.  The  target  height  estimate  is  the  height  that  produces  the  best  match  between  the 
modelled  and  measured  patterns.  Target  elevation  angle  is  derived  from  the  target  height. 
Based  on  the  results  presented,  the  MLE  methods  are  overwhelmingly  superior  to  the  HRDF 
methods.  In  particular,  they  produce  extremely  accurate  estimates  of  target  heights  and 
elevation  angles.  They  are  also  very  tolerant  to  inadequacies  in  array  calibration. 

The  MLE  methods  require  the  development  of  a  practical  propagation  model  that  can 
produce  the  true  target-plus-image  component  of  the  measured  interference  pattern  at  the 
array  elements.  This  task  is  difficult,  because  the  above  component  changes  with  time  and  it  is 
not  always  possible  to  measure  many  key  parameters,  e.g.  temperature  and  moisture  gradients 
at  the  water  surface,  in  battlefield  conditions. 

HRDF  methods  do  not  need  a  propagation  model.  However,  they  require  a  reasonably 
accurately  calibrated  antenna  array  for  implementation  in  LAT  systems.  The  calibration  is 
required  to  compensate  for  unequal  gain  and  phase  responses  at  the  array  elements,  antenna 
tilt,  mutual  coupling  among  the  elements,  the  effect  of  the  ship=s  structure,  and  small  errors  in 
element  positions.  Some  negative  effects  of  calibration  errors  on  direction  estimates  have 
been  discussed  by  Friedlander  [14,15]. 

Several  factors  influence  the  design  of  calibration  methods  for  LAT  arrays.  Firstly,  an 
off-site  calibration  still  requires  an  on-site  refinement.  The  off-site  calibration  cannot 
compensate  for  antenna  tilt  and  the  effect  of  the  ship=s  structure.  Besides,  the  process  of 
moving  the  antenna  to  a  ship  may  change  the  relative  gain  and  phase  responses  at  the  array 
elements.  Secondly,  for  an  on-site  refinement  with  a  near-field  source,  the  effect  of  the  ship=s 
structures  differs  from  what  it  would  be  with  true  targets  many  kilometres  from  the  ships. 
Finite  source  size  and  errors  in  effective  element  positions  may  also  introduce  errors  in  the 
equalisation  of  element  gain  and  phase  responses.  Thirdly,  for  a  refinement  with  a  fixed 
distant  pilot  source,  the  target  and  image  signals  produce  an  interference  pattern  that  changes 
with  time. 
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The  antenna  calibration  methods  in  [5-13]  were  not  designed  for  HRDF  methods.  In 
[8]  and  [10],  for  example,  Lo  and  Litva  only  ensured  accurate  alignment  of  element  positions 
in  the  field.  In  [12]  and  [13],  Bosse,  Turner,  Brookes  and  Riseborough  used  a  calibration 
source  4.6  km  from  the  antenna  array.  They  used  the  statistical  average  of  sea-state  5  data  to 
refine  the  gain  responses.  The  authors  did  not  refine  the  phase  responses. 

An  effective  in-situ  method  to  refine  the  calibration  of  a  LAT  antenna  array  is 
presented  here.  The  method  is  derived  from  a  matrix  calibration  method  developed  by  Hung 
[16,17].  It  requires  a  large  number  of  array  snapshots  generated  by  moving  a  remote  pilot 
source  over  the  range  of  target  elevation  angles  covered  by  the  array.  In  evaluations  with  a 
HRDF  method,  the  calibration  improved  the  accuracy  of  target  elevation  angle  estimates, 
improved  the  resolution  of  target  and  image  elevation  angles,  suppressed  a  bias  in  elevation 
angle  estimates,  and  suppressed  some  anomalous  tracks  between  the  true  target  and  image 
tracks. 

In  the  following  discussions.  Section  2  describes  the  new  calibration  method.  Sections 
3  and  4  evaluate  the  method  with  simulated  and  experimental  data,  respectively,  and  Section  5 
contains  the  conclusions. 
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2.  Calibration  method 


The  calibration  method  requires  reasonably  accurate  elevation  angles  of  the  target  and 
image  signals  in  the  array  snapshots.  A  method  to  calculate  the  image  elevation  angles  will  be 
presented  in  the  analysis  of  experiment  data.  The  following  contains  a  brief  description  of  the 
signal  model  used,  the  assumptions  made,  and  the  computational  procedure. 

2.1  Signal  model 

There  are  N  array  elements.  The  m-th  ideal  array  snapshot  is  given  by 


Xtn  =  aTma(0Tm)  +  ain,a(0|m)  ,  (0 

where  arm  and  aim  are  the  complex  amplitudes  of  the  received  target  and  image  signals, 
respectively,  0Tm  and  0im  are  the  target  and  image  signal  elevation  angles,  and  a(0Tm)  and 
a(0im)  are  the  array  steering  vectors  for  0Tm  and  0|m.  This  expression  assumes  a  perfectly 
calibrated  array,  a  pure  tone  signal  source,  no  interaction  among  the  array  elements,  no  noise, 
and  no  interference  besides  the  image  signal. 

The  corresponding  measured  array  snapshot  is  given  by  ym,  with 
ym  ~  ®Xm  ©m  • 

Matrix  B  has  dimension  NHN  and  depends  on  the  gain  and  phase  responses  at  the  array 
elements,  antenna  tilt,  mutual  coupling  among  the  elements,  the  effect  of  the  ship=s  structure, 
and  errors  in  element  positions.  The  vector  Cm  denotes  the  sum  of  noise  and  unknown 
interfering  signals. 

Let  {yi,y2,AAA,yM}  be  a  set  of  M  array  snapshots.  The  calibration  matrix  calculated 
with  these  snapshots  is  identified  here  as  the  NHN  matrix  C  that  minimises  the  objective 
function 


M  .  M 

S=  3  ♦Cym-Xm*V 

m=l  ni= 


3j*Cym* 


(3) 


subject  to  the  uniqueness  condition 

Cm  =  1.0.  (4) 

The  Xm  is  a  component  of  Cy™  and  the  division  by  the  sum  of  ♦Cym*^  is  designed  to  cancel  the 
effect  of  an  unknown  scaling  factor  introduced  by  C.  Matrix  C  is  proportional  to  B  if  each 
Cm  is  a  null  vector. 
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In  the  convention  used,  >"=  denotes  an  estimate,  superscript  >H=  denotes  a  conjugate- 
transpose,  the  set  {Cpq!  p,q=l,2,^=,N}  denotes  the  coefficients  of  matrix  C,  and  ym,  defined 
as 

y^  =  Cy„,  (5) 


denotes  a  calibrated  snapshot. 


2.2  Conditions 

The  following  conditions  are  assumed  to  be  satisfied: 


1 .  The  array  elements  are  narrowband; 

2.  The  process  that  converts  an  ideal  array  snapshot  to  a  measured  snapshot  can  be 
approximated  by  a  matrix  operation; 

3.  The  calibration  data  are  generated  by  moving  a  remote  pilot  source  over  the  range  of 
elevation  angles  covered  by  the  antenna  array; 

4.  There  are  many  array  snapshots  in  the  calibration  data; 

5.  The  target  elevation  angles  are  given;  and 

6.  A  preliminary  calibration  has  been  carried  out,  so  that  one  can  obtain  reasonably 
accurate  estimates  of  the  image  elevation  angles  from  the  calibration  data  by  curve 
fitting. 

Conditions  3  and  4  are  needed  for  tolerance  to  noise,  unknown  interference,  and 
unknowns  in  the  propagation  medium.  One  can  assume  that  Condition  4  is  satisfied  if,  say, 
every  fifth  snapshot  in  the  calibration  data  is  deleted  and  the  new  calibration  matrix  C=  is 
essentially  the  same  as  C,  so  that  C’*C=  is  approximately  a  constant  times  In. 


2.3  Computation  of  calibration  matrix 

The  calculation  of  calibration  matrix  C  is  based  on  the  observations  presented  below. 

A.  Matrix  C  must  be  calculated  iteratively,  yielding  a  sequence  {Co,Ci,C2,AAA},  where 
Co  is  the  initial  estimate  of  C. 

B.  Let  Ci.i  be  an  estimate  of  C,  {ui„,U2„,}  be  an  orthonormal  basis  constructed  from  the 
array  steering  vectors  a(6Tm)  and  afSim).  and 

Zm  =  Ci.iy„  . 

One  can  calculate  an  estimate  of  Xm  as  the  component  of  in  the  signal  subspace 
spanned  by  a(0Tni)  and  a(0i„,),  i.e. 
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Xm  (U  ImZnOUlm  "^(**2111^)^2111  • 


(7) 


C.  Given  the  NHM  data  matrices  X  and  Y,  where  X=(xi,X2,AAA,Xm)  and 
Y=(yi,y2,AAA,yM)»  one  can  estimate  a  new  C  which  results  in  a  smaller  value  of  the 
objective  function  S  by  solving  the  equation 

CiY  =  X  (8) 

as 

Ci  =  (»^^)(Y^)'  (9) 

and  then  impose  the  uniqueness  condition  given  by  Eq.  (4)  on  Q. 

D.  The  initialisation 

Co  =  In  ,  (10) 

where  In  is  an  NHN  identity  matrix,  is  equivalent  to  using  the  original  data  in  the  first 
iteration,  i.e.  Zni=ym- 

E.  Let 

Di  =  Ci-Ci.,  (11) 

and 

Aj  =  max{*Dpq*:  p,q=l,2,AAA,N}  .  (12) 

One  can  terminate  the  iteration  if  (a)  Ai<A=,  where  A=  is  a  prespecified  small  value, 
or  (b)  i>iniax,  where  imax  is  the  maximum  number  of  iterations  specified. 

F.  If  the  objective  function  has  more  than  one  minimum,  the  matrix  C;  that  gives  the 
smallest  value  of  S  is  the  desired  C. 

The  following  procedure  to  calculate  C  is  based  on  the  above  observations.  There  are 
seven  steps. 

Step  1. 

Assign  imax  and  A=. 

Let  i=l  and  Co=In. 


Step  2. 

For  m=l,2,  AAA,M 
(a)  calculate  Zm  with  (6), 
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(b)  calculate  {uin„U2m}  with  the  Gram-Schmidt  orthonormalization  method 

(c)  calculate  Xm  with  (7). 


Step  3 

Calculate  the  value  of  S  in  the  (i-l)th  iteration  with  (3). 

Step  4. 

Calculate  Ci  with  (9). 

Impose  the  uniqueness  condition  (4)  on  Q. 

Step  5. 

Calculate  A,  with  (11)  and  (12). 

Step  6. 

If  Aj^A—  or  i^imaxj 

Terminate  the  iterations  and  let  C=Ci. 

Else, 

Replace  i  by  i+1  and  return  to  Step  2. 

Step  7. 

Repeat  Steps  2  to  6  with  many  other  initial  values  of  Co.  Study  the  values  of  S  and  the 
coefficients  in  C  as  well  as  the  dependence  of  S  and  Cj  on  Co.  Identify  the  Cj  with  the 
overall  smallest  value  of  S  as  the  desired  C. 


6 


TM  2000-122 


3.  Evaluation  with  simulated  data 


This  evaluation  was  restricted  to  simulated  data  with  unequal  element  gain  and  phase 
responses  only.  It  was  designed  to  show  the  adverse  effects  of  unequal  responses  on  elevation 
angle  estimates  calculated  with  a  HRDF  method  and  illustrate  how  calibration  could  alleviate 
these  effects.  The  radar,  target,  and  propagation  parameters  were  chosen  to  produce  results 
similar  to  those  of  the  experimental  evaluation  in  the  next  section. 

3.1  Data  Generation 

The  receive  anteima  was  a  vertical,  uniform,  linear  array  with  N=8  elements  spaced  at 
a  distance  of  12.5  cm.  It  was  centred  at  a  height  hR=6.7  m  above  water.  The  operating 
frequency  was  9.6  GHz,  yielding  a  signal  wavelength  of  X,=0.03125  m  and  an  element  spacing 
in  wavelengths  of  d=4.0.  The  principal  range  of  elevation  angles  was  -7.2E  to  7.2E. 

The  ideal  snapshot  Xm  was  generated  using  a  flat  earth  model,  Eq.  (1),  and 


©Tm  =  tan’'((hTm-hR)/RTni)  , 

(13) 

Gin,  =  tan''((-hTm-hR)/RTni)  , 

(14) 

arm  =  exp(j<|)Tm)  , 

(15) 

aim  =  paTmexp(j54,m) , 

(16) 

“  *2llARn/A.  , 

(17) 

ARm  “  2hxmhR/RTiTi  > 

(18) 

a(0)  =  (l,exp[j2ndsin0],=S5,exp[j27td(N-l)sin0])^ . 

(19) 

The  parameter  hxm  was  the  target  height,  Rim  was  the  target  range,  SNRt  was  the  target  SNR 
at  the  elements,  (|)Tm  was  a  random  phase,  p  was  the  reflection  coefficient  at  the  water  surface, 
Arto  was  the  transit  path  difference  between  the  target  and  image  signals,  and  was  the 
phase  correction  for  this  difference.  The  p  was  assigned  as 

p  =  -0.9 ,  (20) 

corresponding  to  0.92  dB  attenuation  and  180E  phase  change.  In  (19),  the  first  element  was 
the  lowest  element  in  the  array  and  0  represented  Gm,  or  Oim- 

The  measured  snapshot  ym  was  generated  with  (2)  and  either 

B  =  diag{gi,gi,AAA,gN}  (21a) 


or 


TM  2000-122 


7 


(21b) 


B  =  diag{gi,g2,AAA,gN}  . 

The  gn's  were  complex  amplitudes  of  element  gain  and  phase  responses  and  had  the  values 
given  in  Table  1.  These  amplitudes  were  chosen  to  be  the  same  as  those  calculated  from  the 
experimental  data  in  the  next  section.  The  components  of  Cm  in  ym  were  generated  as 
uncorrelated  complex  Gaussian  variates  with  zero  mean  and  unit  variance. 


n 

Gain*g„* 

Phase  (°) 

1 

1.000 

0.0 

2 

1.044 

-14.0 

3 

0.988 

-7.8 

4 

1.098 

-16.6 

5 

0.978 

-3.8 

6 

1.047 

1.4 

7 

0.978 

8.2 

8 

0.966 

13.9 

Tablel.  Relative  element  gain  jgn!  and  phase  ang(gn)  responses  in  simulation. 


Three  files  were  produced.  File  1  used  the  B  given  by  (21a)  and  had  M=1971 
snapshots.  The  target  was  a  transmitter  with  SNRt  =30  dB  at  a  fixed  range  RTn,=8.0  km.  It 
moved  up  from  height  hTi=7.0  m  to  hTM=992.0  m  (0ti=O.OOE,  6ii=-0.01E,  0tm=7.O2E,  and 
0,m=-7.16E)  while  the  snapshots  were  generated  at  intervals  of  0.5  m.  File  2  was  generated 
with  the  same  B  and  had  M=451  snapshots.  The  target  had  SNTt=20  dB  and  a  fixed  height 
hTm=46.0  m.  It  approached  the  radar  from  range  Rti=20.0  km  to  Rtm=2.0  km  while  the 
snapshots  were  generated  at  range  intervals  of  4.0  m.  File  3  was  generated  with  the  B  given 
by  (21b).  The  other  parameters  were  the  same  as  those  in  File  2. 

3.2  Calibration  matrix 

File  1  was  used  to  calculate  C.  In  Step  1,  the  assignments  were  A==0.00001  and 
iina!c=400.  In  Step  7,  the  C  calculated  with  Co=In  was  treated  as  the  output  of  a  coarse  search. 
It  was  refined  by  repeating  Steps  2  and  6  fifteen  times.  In  each  repetition,  the  new  Co  was 
obtained  by  adding  complex  amplitudes  with  zero  mean  and  0.04  standard  deviation  to  the  C 
obtained  in  the  coarse  search.  The  step  also  included  many  other  coarse  searches  followed  by 
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refinements.  In  each  coarse  search,  Co  was  obtained  by  adding  complex  amplitudes  with  zero 
mean  and  0.3  standard  deviation  to  every  coefficient  of  In. 

The  C  matrices  at  the  end  of  these  coarse-plus-refinement  searches  were  essentially 
the  same  and  were  approximately  diagonal.  All  those  selected  for  further  evaluation  produced 
similar  results.  The  values  of  S  for  these  matrices  were  between  0.000407  and  0.000445.  In 
the  Cj  with  S=0 .000407,  all  off-diagonal  coefficients  did  not  exceed  0.001  in  magnitude  and 
that  the  product  CjB  was  an  8H8  identity  matrix  when  only  the  first  three  decimal  places  were 
retained. 

In  the  calculation  with  Co=In,  the  Aj  in  Step  5  initially  decreased  monotonically  from 
0.2634584  at  i=l  to  0.0000109  at  i=181.  Then,  it  fluctuated  with  amplitudes  less  than 
0.0000001  and  decreased  further  to  less  than  0.0000100  at  i=198.  The  corresponding  value  of 
S  decreased  monotonically  from  0.025769  at  i=l  to  0.000435  at  i=198.  Refinement  produced 
a  new  C  with  S=0.000434. 

3.3  Elevation  angle  estimates 

The  TKSVD  method  [2,3]  was  used  to  calculate  six  candidate  direction  estimates 
from  each  snapshot  in  File  2,  using  a  sixth  order  linear  prediction  filter  and  assuming  that  two 
signal  sources  were  present.  The  two  candidate  direction  estimates  with  the  largest  powers 
were  identified  as  target  and  image  elevation  angle  estimates.  Fig.  1  shows  the  estimates, 
plotted  as  points,  before  (B,  upper  half)  and  after  (A,  lower)  calibration.  Also  plotted  are  the 
target  (solid  line)  and  image  (dashed  line)  tracks. 

There  are  several  notable  observations.  Firstly,  the  TKSVD  estimates  form  two 
visually  identifiable  tracks:  an  upper  track  for  the  target,  and  a  lower  track  for  the  image. 
Secondly,  there  is  a  bias  in  the  TKSVD  estimates  before  calibration,  because  the  tracks 
produced  by  the  TKSVD  estimates  are  above  the  true  tracks  at  ranges  less  than  5.0  km.  This 
bias  is  not  observable  after  calibration.  Thirdly,  before  calibration,  there  are  estimates  that 
produce  anomalous  tracks  between  the  true  target  and  image  tracks  at  ranges  4.4,  5.6,  7.9  and 
13.2  km.  After  calibration,  the  anomalous  tracks  are  replaced  by  horizontal  tracks  with 
approximately  zero  elevations.  Finally,  outside  the  4.4,  5.6,  7.9  and  13.2  km  regions,  the 
individual  estimates  are  closer  to  the  true  target  and  image  tracks  after  calibration. 

Fig.  2  was  calculated  with  File  3,  which  was  generated  with  the  matrix  B  given  by 
(21b).  Only  the  elevation  angle  estimates  before  calibration  are  presented,  because  G  could 
not  be  used  to  calibrate  this  file.  Compared  with  the  upper  half  of  Fig.  1,  both  the  bias  in 
elevation  angle  estimates  and  the  gradients  of  the  anomalous  tracks  have  changed  signs.  These 
sign  changes  show  that  equalisation  of  element  phase  responses  is  very  important,  because 
these  responses  play  a  dominant  role  in  the  generation  of  bias  and  anomalous  tracks. 

A  study  of  the  film's  in  Files  2  and  3  showed  that  the  anomalous  tracks  were  centred  at 
ranges  where  the  target  and  image  signals  had  approximately  the  same  phase  at  the  centre  of 
the  array.  After  calibration,  each  snapshot  close  to  these  ranges  produced  one  elevation  angle 
estimate  with  G.O.OE  plus  one  estimate  with  e.2.7E  or  -2.7E.  The  0.V2.7E  estimates 
correspond  to  noise  pole  positions  in  the  TKSVD  method. 


TM  2000-122 


(deg)  i  (deg) 


Range  (km) 


Range  (km) 


Figure  1.  TKSVD  elevation  angle  estimates  (points),  target  track  (solid  iine),  and  image  track  (dashed 
line)  obtained  from  simulated  data  File  2  before  (B,  upper  half)  and  after  (A,  lower  half)  calibration. 
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(deg) 


Range  (km) 

Figure  2.  TKSVD  elevation  angle  estimates  (points),  target  track  (solid  line),  and  image  track  (dashed 
line)  obtained  from  simulated  data  File  3  before  (B)  calibration. 
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4.  Evaluation  with  experimental  data 

The  experimental  data  sets  in  this  evaluation  were  generated  for  the  evaluation  of 
LAT  methods  only.  They  were  used  to  evaluate  the  calibration  method  because  more  suitable 
data  were  unavailable. 

4.1  Experimental  data 

The  experimental  data  were  measured  with  an  X-band  LAT  radar  in  1981  by  Standard 
Telecommunication  Laboratories  (STL)  at  the  Fraser  Gunnery  Range,  near  Portsmouth, 
United  Kingdom.  This  radar  had  separate  vertically  polarised  transmit  and  receive  antenna 
arrays  mounted  side  by  side  on  a  pedestal  that  could  be  rotated  about  a  vertical  axis.  Between 
them  was  an  optical  tracker  to  provide  an  independent  measure  of  the  target  elevations.  A 
detailed  description  of  the  experiment  was  given  by  Pearson  and  Waddoup  [18]  and  a  brief 
description  is  given  in  [4].  Here,  it  is  sufficient  to  know  that,  like  the  set-up  in  the  simulation 
evaluation,  the  transmit  antenna  operated  at  9.6  GHz,  the  receive  antenna  was  a  vertical, 
uniform,  linear  array  with  N=8  elements  spaced  at  a  distance  of  d=4.0  wavelengths,  and  the 
principal  range  of  target  and  image  elevation  angles  is  -7.2E  to  7.2E. 

This  author  had  access  to  five  data  files.  Each  file  was  generated  with  a  Hunter  or 
Canbera  aircraft  flying  towards  the  radar.  It  contained  snapshots  of  the  array  element  outputs 
as  well  as  the  aircraft  ranges  and  elevations.  Table  2  gives  the  file  identifications,  dates  of 
measurement,  aircraft  types,  number  of  array  snapshots  in  the  file,  and  the  receive  antenna 
heights.  STL  File  3  originally  had  440  snapshots.  Only  the  first  265  were  retained,  because  the 
close-range  snapshots  had  large  range  and  elevation  errors.  Similarly,  STL  File  5  originally 
had  510  snapshots  and  only  the  first  405  were  retained. 


STL  File 

1 

2 

3 

4 

5 

ID 

250 

252 

311 

312 

313 

1981  Date 

June  18 

June  18 

June  25 

June  25 

June  25 

Aircraft 

Hunter 

Hunter 

Canberra 

Canberra 

Canberra 

Snapshots 

330 

320 

265 

415 

405 

Nominal  hr  (m) 

46 

46 

61 

61 

61 

Nominal  ha  (m) 

6.7 

6.7 

mm 

7.9 

Bias  (deg) 

0.23 

0.23 

0.18 

0.18 

0.18 

Table  2.  STL  low-angle  tracking  data  files 
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4.2  Elevation  estimates  with  uncalibrated  data 

The  TKSVD  method  was  used  to  calculate  two  elevation  angle  estimates  from  each 
snapshot  in  the  files.  The  results  obtained  with  STL  Files  1  to  5  are  plotted  in  the  upper  halves  (B) 
of  Figs.  3  to  7,  respectively.  Also  plotted  are  the  target  tracks  given  by  the  optical  tracker  (solid 
lines).  There  are  two  notable  observations.  Firstly,  the  results  look  similar  to  those  in  the  upper 
half  of  Fig.  1.  Secondly,  the  bias  in  the  TKSVD  estimates  depends  on  the  date  of  the  experiment. 

The  first  observation  is  expected,  because  the  parameters  in  the  simulated  data  were 
chosen  to  produce  results  similar  to  the  STL  results.  The  second  observation  was  first  noted  in 
[18],  where  it  was  stated  that  the  exact  reason  was  unknown. 

The  bias  in  a  given  file  is  denoted  by  0b  and  is  given  in  Table  2.  It  is  defined  as  the 
up-shift  that  must  be  applied  to  the  optical  track  to  make  it  fit  the  TKSVD  target  elevation 
angle  estimates.  The  value  is  0b=O.23E  for  Files  1  and  2  generated  on  June  18  and  0b=O.18E 
for  Files  3  to  5  generated  on  25. 

4.3  Calibration  matrix 

Files  1  and  2  had  the  same  bias  and  their  optical  tracks  fitted  the  bias-removed 
TKSVD  target  elevation  estimates  very  well.  They  were  merged  to  produce  a  new  file  with 
650  snapshots.  This  new  file  was  used  to  calculate  C,  although  Conditions  3  and  4  in  Section 
2.2  were  not  satisfied.  In  Step  1,  the  assignments  were  A==0.0001  and  imax=300-  In  Step  2,  the 
target  and  image  elevation  angles  were  identified  as  the  elevation  angles  in  the  target  and 
image  tracks,  respectively.  In  Step  3,  was  highly  ill  conditioned,  because  the  target 
elevations  did  not  exceed  1.5®  and  was  well  below  the  7.2°  upper  limit  of  principal  target 
elevation  angles.  Consequently,  the  (Y^)  '  in  Eq.  (9)  was  replaced  by  the  Moore-Penrose 
pseudo-inverse  calculated  with  the  singular  value  decomposition  method.  Only  the  four 
largest  singular  values  were  not  equated  to  zeroes,  because  this  number  produced  the  best 
results  after  calibration.  Step  7  was  omitted.  In  the  result,  Rank(C)=4  and  C  was  not 
symmetric,  not  diagonal,  and  not  even  approximately  diagonal. 

The  intermediate  results  showed  that  the  Aj  in  Step  5  initially  decreased 
monotonically  from  0.73656  for  i=l  to  0.00202  for  i=63.  Then,  it  fluctuated  with  amplitudes 
less  than  0.00003  and  decreased  further  to  0.00069  for  i=300.  The  values  of  S  was  0.03222 
for  i=l  and  0.00351  for  i=2.  It  then  increased  slowly  between  i=2  and  i=53  to  0.00462  before 
decreasing  monotonically  to  0.00436  for  i=300.  Between  iterations  i=2  and  i=53,  C  gradually 
acquired  the  properties  noted  in  Section  4.4  below.  The  decrease  between  i=53  and  i=300 
corresponded  to  a  refinement  of  C. 

4.4  Elevation  estimates  with  calibrated  data 

The  STL  files  were  calibrated  with  C,  using  Eq.  (5).  The  new  TKSVD  elevation  angle 
estimates  are  plotted  in  the  lower  halves  of  Figs.  3  to  7.  Evidently,  the  calibration  has 
improved  the  accuracy  of  the  estimates,  improved  the  resolution  of  target  and  image  elevation 
angles,  suppressed  the  anomalous  tracks  between  the  target  and  image  tracks,  suppressed  the 
bias  in  STL  Files  1  and  2,  and  reduced  the  bias  in  Files  3  to  5. 


TM  2000-122 


13 


The  degree  of  improvement  was  better  than  expected,  because  the  number  of 
snapshots  in  the  merged  file  was  rather  small,  the  SNRt=s  in  snapshots  were  significantly 
lower  than  what  one  would  get  with  a  pilot  source,  and  the  largest  target  elevation  angle  was 
significantly  below  7.2®.  The  bias  in  Files  3  to  5  were  only  reduced,  because  the  bias  in  these 
files  was  different  from  the  bias  in  the  data  used  to  calculate  C. 

4.5  Relative  element  gains  and  phases 

The  vector  g(0),  given  by 

g(0)  =  cVe), 

provides  a  measure  of  the  relative  gain  and  phase  responses  at  the  array  elements  for  signal 
elevation  0.  Table  I  gives  the  responses,  relative  to  the  first  element,  for  0=O.OE.  These  values 
were  used  in  the  generation  of  the  simulated  data  used  in  Section  3. 


(22) 
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(deg)  a  (deg) 


Range  (km) 


Figure  3.  TKSVD  elevation  angle  estimates  (points),  optical  target  track  (solid  line),  and  image  track 
(dashed  line)  obtained  from  STL  File  1  before  (B)  and  after  (A)  calibration. 
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Figure  4.  TKSVD  elevation  angle  estimates  (points),  optical  target  track  (solid  line),  and  image  track 
(dashed  line)  obtained  from  STL  File  2  before  (B)  and  after  (A)  calibration. 
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Figure  5.  TKSVD  elevation  angle  estimates  (points)  and  optical  target  track  (solid  line)  obtained  from 

STL  File  3  before  (B)  and  after  (A)  calibration. 
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Figure  6.  TKSVD  elevation  angle  estimates  (points)  and  optical  target  track  (solid  line)  obtained  from 

STL  File  4  before  (B)  and  after  (A)  calibration. 
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Figure  7.TKSVD  elevation  angle  estimates  (points)  and  optical  target  track  (solid  line)  obtained  from  STL 

Fiie  5  before  (B)  and  after  (A)  calibration. 
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5.  Conclusion 


This  author  has  described  an  in-situ  method  to  calibrate  an  X-band  antenna  array  for 
tracking  low-elevation  targets  with  HRDF  methods.  The  simulation  and  experimental 
evaluations  showed  that  the  calibration  (a)  improved  the  accuracy  of  target  elevation  angle 
estimates  calculated  with  the  TKSVD  method,  (b)  improved  the  resolution  of  target  and  image 
elevation  angles,  (c)  suppressed  a  bias  in  elevation  angle  estimates,  and  (d)  suppressed  the 
anomalous  tracks  between  the  true  target  and  image  tracks. 
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